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ABSTRACT 

We present a sample of low-resolution 5-38 /im Spitzer IRS spectra of the inner few square kiloparsecs 
of 59 nearby galaxies spanning a large range of star formation properties. A robust method for 
decomposing mid-infrared galaxy spectra is described, and used to explore the behavior of PAH 
emission and the prevalence of silicate dust extinction. Evidence for silicate extinction is found in 
~l/8 of the sample, at strengths which indicate most normal galaxies undergo Ay <3 magnitudes 
averaged over their centers. The contribution of PAH emission to the total infrared power is found 
to peak near 10% and extend up to ^20%, and is suppressed at metallicities Z < z7 /4, as well as 
in low-luminosity AGN environments. Strong inter-band PAH feature strength variations (2-5 x) are 
observed, with the presence of a weak AGN and, to a lesser degree, increasing metallicity shifting power 
to the longer wavelength bands. A peculiar PAH emission spectrum with markedly diminished 5-8 fiui 
features arises among the sample solely in systems with relatively hard radiation fields harboring 
low-luminosity AGN. The AGN may modify the emitting grain distribution and provide the direct 
excitation source of the unusual PAH emission, which cautions against using absolute PAH strength to 
estimate star formation rates in systems harboring active nuclei. Alternatively, the low star formation 
intensity often associated with weak AGN may affect the spectrum. The effect of variations in the 
mid-infrared spectrum on broadband infrared surveys is modeled, and points to more than a factor 
of two uncertainty in results which assume a fixed PAH emission spectrum, for redshifts z = — 2.5. 
Subject headings: galaxies: active — galaxies: ISM — infrared: galaxies — techniques: spectroscopic 



1. INTRODUCTION 

The mid-infrared (MIR) spectrum of star-forming 
galaxies is dominated by strong emission features gen- 
erally attributed to Polycyclic Aromatic Hydrocarbons 
(PAHs). PAHs are aromatic molecules ubiquitous in the 
interstellar medium (ISM) of our own Galaxy, and found 
in most nearby galaxies with ongoing or recent star for- 
mation. These large molecules, typically a few A in di- 
ameter and harboring up to a few hundred carbon atoms, 
are now commonly believed to be responsible for the 
prominent suite of broad emission bands in the MIR be- 
tween 3-19 /im (Leger & Puget 1984; Allamandola et al. 
1985). PAH molecules efficiently absorb ultraviolet and 
optical photons, and in many environments supply the 
bulk of the photo-electrons responsible for heating inter- 
stellar gas (e.g. Weingartner & Draine 2001). 

The mid-infrared PAH bands can be very strong, both 
compared to their underlying dust grain continuum, and 
even to the reprocessed bolometric emission of larger 
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dust grains. In galaxies with intense star formation, 
up to 20% of the total infrared luminosity is emitted 
in the PAH bands alone, and, longward of the Lyman 
and Balmer breaks, they are among the most prominent 
spectral features with sufficient contrast for photometric 
redshift determination. In deep infrared surveys, they 
can induce peaks in the flux-limited number density of 
15-24 /mi-detected galaxies at redshifts z >1, as they 
move through fixed filter passbands (e.g. Barvainis et al. 
1999; Le Floc'h et al. 2005; Caputi et al. 2006). Thus 
they are of fundamental importance in modeling deep in- 
frared survey results, and especially in constraining the 
role of dust-obscured star formation at early epochs. Yet 
the behavior or even presence of PAHs in the distant 
universe is largely unexplored. Despite tantalizing indi- 
cations of PAH emission at redshift beyond z>2 (Yan 
et al. 2005; Lutz et al. 2005; Huang et al. 2006), other 
surveys probing the mid-infrared properties of optically 
faint distant galaxy samples have uncovered mid-infrared 
spectra dominated by continuum emission, cut by deep 
silicate absorption (Houck et al. 2005). 

A prevailing difficulty in studying PAH emission at 
high redshifts is the relative ambiguity of the 7-13 /im 
spectrum of star-forming galaxies, with the broad PAH 
features at 7.7 /im and 11.3/zm almost perfectly brack- 
eting the 10 /im silicate absorption trough. Silicate ab- 
sorption can be pronounced in some galaxies' spectra, 
including high redshift objects dominated by AGN emis- 
sion. At low signal, or in broad photometric bands, dis- 
criminating between these two physically distinct but ob- 
servationally similar signatures can be challenging. In- 
terpretation of these spectral signatures, and of broad- 
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band properties of even fainter and more distant systems, 
has relied predominantly on a small set of local template 
spectra developed primarily from Infrared Space Obser- 
vatory (ISO) observations of bright nearby systems. 

The Spitzer Infrared Nearby Galaxies Survey (SINGS, 
Kennicutt et al. 2003) provides an exceptional sample 
in which to explore the global and spatially resolved be- 
havior of these fundamental emission bands over a much 
broader range of physical environments than has been 
possible before. We focus here on the global behavior of 
the PAH emission spectrum and its underlying contin- 
uum. A future paper will focus on the spatially resolved 
behavior of PAH emission in the centers and extended 
disks of galaxies in the SINGS sample. 

2. BACKGROUND ON THE EMISSION PROCESSES 

The "unidentified infrared bands" were first discovered 
in 1973 in the spectra of Galactic planetary nebulae by 
Gillett ct al. PAHs were later proposed as the origin of 
these features (Leger & Puget 1984; Allamandola et al. 
1985), and invoked to explain the 12 /im infrared excess 
seen by IRAS in the diffuse cirrus of the Milky way 
(Puget et al. 1985). Although the carriers of the mid- 
infrared emission bands are now commonly attributed 
to the broad class of PAH molecules, other interpreta- 
tions remain consistent with at least some aspects of the 
bands' strength, shape, and location (see Li & Drainc 
2001, and references therein). Although the only cer- 
tainty is that they arise primarily in material with sp 2 
hybridized carbon-carbon bonds by infrared fluorescence, 
we adopt the convention of referring to them as PAH 
emission bands; it should be kept in mind that the emit- 
ters could encompass a larger group of individual species. 
In this same spirit, we refer to PAHs interchangeably as 
both molecules and small grains. 

The ISO Key Project on Nearby Galaxies (Helou et al. 
2000) explored the 2-11 /zm spectrum of a sample of nor- 
mal galaxies. ISO was limited to relatively high sur- 
face brightness sources, whereas SINGS extends to much 
lower luminosity and surface brightness, and is exploring 
a much broader range of local interstellar environment. 
A surprising result of the Key Project was the relative 
uniformity of the PAH emission spectrum over a range 
of different galaxy environments. Lu et al. (2003) found 
band to band variations of at most 15% out to 11 /im. 
A simple physical argument which can be invoked to ex- 
plain this apparent uniformity relates to the underlying 
emission mechanism of small grains. Since PAH grains 
are transiently excited by single photons to high temper- 
atures (higher than their equilibrium temperatures; see 
Leger & Puget 1984; Draine & Anderson 1985; Allaman- 
dola et al. 1985; Draine & Li 2001), the emission spec- 
trum produced by a single grain population is relatively 
insensitive to the strength and hardness of the exciting 
radiation field (Leger & Puget 1984), and the resulting 
PAH spectrum might therefore be expected to be quite 
uniform. However, a variety of mechanisms are known 
to affect to various degrees the strength of the individual 
PAH bands, including ionization state, distribution of 
grain sizes, de-hydrogenation, and multi-photon heating 
in intense radiation fields (see Li 2004, for a review). The 
possible destruction mechanisms of PAH grains are var- 
ied as well, including photo-thermo (Leger et al. 1989; 
Guhathakurta & Draine 1989) and Coulomb dissocia- 



tion, sputtering, and dissociation in shocks. 

Evidence in our own Galaxy (e.g. Joblin et al. 1996; 
Verstraete et al. 1996; Hony et al. 2001; Peeters et al. 
2004; Bregman & Temi 2005) and resolved regions in 
nearby galaxies like the LMC and SMC (Contursi et al. 
2000; Vermeij et al. 2002), indicate that the PAH emis- 
sion spectrum is variable on small scales, and can depend 
sensitively on the local gas density and details of the inci- 
dent radiation field. Variations can occur in the relative 
strengths, widths, and wavelength centers of individual 
emission bands (e.g. Peeters et al. 2002). 

It is not yet known whether the primary processes reg- 
ulating PAH emission are chemical, mechanical, or radia- 
tive, or whether they are expressed through inhibited or 
enhanced grain formation rates, selective destruction in 
differing environments, or varying efficiencies of photo- 
excitation. Given the predicted and observed diversity 
in the form of the PAH emission spectrum within indi- 
vidual objects in the Galaxy, the question which remains, 
and which will be explored here, is whether the natural 
variations found on small scales vanish when averaged 
over real distributions of emitting environments in the 
disks and nuclei of galaxies. 

3. OBSERVATIONS AND REDUCTION 

The SINGS survey includes uniform 5-38 /im spectral 
mapping with the Infrared Spectrograph (IRS, Houck 
et al. 2004). Full observational details regarding the 
SINGS spectral mapping program can be found in Ken- 
nicutt et al. (2003). Briefly, a low- resolution (R ~ 50 — 
100), 5-15 /xm spectral map of approximate size 55"x34" 
is centered on the nucleus of each target. Larger, 15- 
38 /im radial strip maps extending to roughly half the 
optical radius and covering 50"x3.5'-15.4' are also ob- 
tained. The total effective integration time per pixel was 
uniform, at 28 seconds in the short wavelength module, 
and 2 minutes in the long wavelength module. The IRS 
spectra were supplemented for all sample targets with 
scan-map images at 24, 70, and 160 /im obtained with 
the MIPS instrument (Rieke et al. 2004). 

3.1. Sample 

A subset of 59 galaxies from the SINGS sample was 
selected, and their relevant parameters are given in Ta- 
ble 1. The sub-sample spans a range in luminosity and 
infrared/optical ratio similar to the complete sample, 
and includes galaxies from all represented morphologi- 
cal types. The only requirements for inclusion in the 
sub-sample were full 5-38 /im spectral coverage, and suf- 
ficient signal-to-noise to put useful upper limits on PAH 
emission (the continuum must be detectable). Among 
the galaxies not included were an individual source with 
anomalous silicate extinction (Roussel et al. 2006), a 
source with strong silicate emission, and several unde- 
tected low-luminosity dwarf galaxies. 

3.2. MIPS Data Reduction and Matched Fluxes 

All MIPS images were reduced as described in Dale 
et al. (2005), and global MIPS fluxes were taken from 
that work, modified slightly (<10%) to account for the 
most recent MIPS calibration factors. Obtaining total 
infrared luminosities (Ltir) matched to the apertures 
listed in Table 1, which are in most cases smaller than 
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TABLE 1 
Sample Parameters 
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0.3 


0.06 


3.22 
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41.9 
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1.88 


NGC0024 


Hn 


7.3 
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35.1 x 29.6 


1.32 


-39.5 


0.9 


0.17 


0.76 


NGC0337 


Hn 


22.4 


00:59:50.20 


-07:34:45.8 


42.6 x 25.9 


12.99 


149.0 


60.4 


0.37 


3.40 


NGC0584 


LINER 


20.0 


01:31:20.90 
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18.5 x 18.5 


3.20 


-32.5 


0.7 


0.12 


0.03 


NGC0628 


Hn 


7.3 


01:36:41.60 


+15:47:00.0 


53.6 x 29.6 


1.99 


151.8 


2.7 


0.03 


0.90 


NGC0855 


Hn 


9.6 


02:14:03.70 


+27:52:38.4 


18.5 x 24.1 


0.95 


156.1 


2.1 


0.55 




NGC0925 


Hn 


9.1 


02:27:17.25 


+33:34:41.6 


27.8 x 29.6 


1.61 


156.2 
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0.05 


1.17 


NGC1097 


LINER 


17.1 


02:46:18.86 


-30:16:27.2 


37.0 x 31.4 


7.99 


146.7 


249.3 


0.38 


6.19 


NGC1266 


LINER 


30.0 


03:16:00.71 


-02:25:36.9 


22.2 x 20.4 


9.58 


157.7 


141.9 


0.54 


16.24 


NGC1291 


LINER 


10.8 


03:17:18.59 


-41:06:28.0 


38.8 x 27.7 


2.97 


-78.5 


4.2 


0.11 


0.17 


NGC1316 


LINER 


24.3 


03:22:41.68 


-37:12:29.4 


22.2 x 25.9 


7.99 


149.0 


27.7 


0.28 


0.27 


NGC1404 


LINER 


18.5 


03:38:51.95 


-35:35:39.1 


25.9 x 22.2 


4.60 


-79.7 


0.6 


0.18 


0.02 


NGC1482 


Hn 


23.2 


03:54:38.88 


-20:30:07.1 


27.7x29.6 


10.40 


-171.4 


299.2 


0.68 


38.04 


NGC1512 


LINER 


11.8 


04:03:54.17 


-43:20:54.4 


22.2 x 25.9 


1.89 


141.9 


8.7 


0.21 


1.62 


NGC1566 


Seyfert 


20.3 


04:20:00.33 


-54:56:16.6 


25.9 x 24.1 


6.06 


129.0 


47.2 


0.08 


2.06 


NGC1705 


Hn 


5.1 


04:54:14.50 


-53:21:36.4 


14.8 x 14.8 


0.13 


120.9 


0.1 


0.19 


0.41 


NGC2403 


Hn 


3.2 


07:36:49.95 


+65:36:03.5 


51.8 x 33.3 


0.42 


30.4 


0.8 


0.02 


0.62 


NGC2798 


Hn 


26.2 


09:17:22.80 


+41:59:59.4 


25.9 x 22.2 


9.26 


23.2 


246.4 


0.66 


6.69 


NGC2841 


Seyfert 


14.1 


09:22:02.50 


+50:58:34.1 


22.2 x 22.2 


2.30 


38.6 


4.0 


0.03 


0.11 


NGC2915 


Hn 


3.8 


09:26:10.03 


-76:37:32.2 


18.5 x 18.5 


0.11 


8.4 


0.1 


0.33 




NGC2976 


Hn 


3.6 


09:47:15.22 


+67:55:00.3 


38.8 x 27.8 


0.32 


32.5 


0.7 


0.08 


0.93 


NGC3049 


Hn 


23.9 


09:54:49.59 


+09:16:18.1 


37.0 x 20.4 


10.11 


15.1 


32.7 


0.61 


2.18 


NGC3184 


Hn 


11.1 


10:18:16.90 


+41:25:24.7 


14.8 x 13.0 


0.56 


33.9 


2.8 


0.03 


1.28 


NGC3190 


LINER 


20.9 


10:18:05.63 


+21:49:54.2 


27.7 x 31.5 


8.96 


18.5 


30.9 


0.37 


0.58 


NGC3198 


LINER 


13.7 


10:19:54.84 


+45:32:58.7 


22.2 x 22.2 


2.17 


36.0 


15.6 


0.17 


5.46 


NGC3265 


Hn 


23.2 


10:31:06.80 


+28:47:45.6 


27.8 x 22.2 


7.76 


22.1 


25.2 


0.71 


1.81 


NGC3351 


Hn 


9.3 


10:43:57.72 


+11:42:13.5 


40.7 x 29.6 


2.46 


-162.3 


28.5 


0.36 


3.20 


NGC3521 


LINER 


10.1 


11:05:48.58 


-00:02:07.3 


53.7 x 31.4 


4.01 


-162.1 


28.3 


0.10 


1.21 


NGC3621 


LINER 


6.6 


11:18:16.51 


-32:48:49.3 


53.6 x 31.4 


1.75 


-145.9 


6.1 


0.07 


1.71 


NGC3627 


Seyfert 


9.4 


11:20:15.04 


+12:59:29.0 


44.4 x 31.5 


2.89 


22.8 


23.5 


0.09 


0.82 


NGC3773 


Hn 


11.9 


11:38:12.98 


+12:06:45.8 


24.1 x 24.0 


1.92 


-159.6 


4.1 


0.53 


0.82 


NGC3938 


LINER 


13.3 


11:52:49.32 


+44:07:13.6 


49.9 x 29.6 


6.19 


36.6 


11.2 


0.10 


0.75 


NGC4125 


LINER 


22.9 


12:08:05.84 


+65:10:29.5 


24.0 x 24.1 


7.13 


173.4 


5.1 


0.32 


0.06 


NGC4254 


Hn 


16.6 


12:18:49.57 


+14:24:57.5 


40.7 x 24.1 


6.33 


18.9 


56.7 


0.11 


1.59 


NGC4321 


LINER 


14.3 


12:22:54.87 


+15:49:19.2 


31.4 x 25.9 


3.92 


24.6 


54.6 


0.16 


3.39 


NGC4450 


LINER 


16.6 


12:28:29.71 


+17:05:08.7 


48.1 x 29.6 


9.20 


31.4 


9.8 


0.20 


0.23 


NGC4536 


Hn 


14.4 


12:34:27.03 


+02:11:16.5 


25.9 x 24.0 


3.06 


21.4 


96.3 


0.73 


9.12 


NGC4552 


Seyfert 


15.9 


12:35:39.88 


+12:33:23.3 


25.9 x 24.0 


3.71 


-160.9 


0.7 


0.25 


0.01 


NGC4559 


Hn 


10.3 


12:35:57.58 


+27:57:34.2 


51.8 x 33.3 


4.29 


40.9 


8.8 


0.13 


0.79 


NGC4569 


Seyfert 


16.6 


12:36:49.76 


+13:09:45.5 


25.9 x 22.2 


3.72 


25.4 


45.3 


0.31 




NGC4579 
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16.6 


12:37:43.53 


+11:49:03.8 


22.2 x 18.5 


2.65 


19.6 


17.1 


0.12 


0.45 
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Seyfert 


9.3 


12:39:59.56 
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0.08 
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Hn 
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37.0 x 29.6 


2.16 


46.3 


2.4 


0.39 


0.68 


NGC4631 


Hn 


8.1 


12:42:07.80 


+32:32:34.6 


50.0 x 33.3 


2.54 


-175.9 


32.6 


0.12 


5.21 


NGC4725 


Seyfert 


11.9 


12:50:26.59 


+25:30:01.2 


31.4 x 20.4 


2.13 


31.2 


3.6 


0.04 


0.28 


NGC4736 


LINER 


5.0 


12:50:53.15 


+41:07:14.4 


53.7 x 33.3 


1.05 


-161.9 


12.2 


0.18 


0.35 


NGC4826 


Seyfert 


5.0 


12:56:43.59 


+21:40:58.0 


29.6 x 31.5 


0.54 


30.2 


12.3 


0.34 


2.38 


NGC5033 


Seyfert 


14.8 


13:13:27.32 


+36:35:35.2 


49.9 x 31.5 


8.11 


53.1 


61.4 


0.24 


1.95 


NGC5055 


LINER 


7.8 


13:15:49.35 


+42:01:45.7 


38.8 x 29.6 


1.65 


22.9 


16.9 


0.08 




NGC5194 


Seyfert 


7.8 


13:29:52.80 


+47:11:43.5 


55.5 x 31.5 


2.48 


61.9 


32.8 


0.08 


2.36 


NGC5195 


Seyfert 


8.0 


13:29:59.50 


+47:15:56.7 


33.3 x 25.9 


1.30 


24.3 


16.4 


0.75 


2.66 


NGC5713 


Hn 


29.4 


14:40:11.38 


-00:17:24.2 


44.4 x 29.6 


26.72 


18.7 


318.9 


0.47 


5.99 


NGC5866 


LINER 


15.1 


15:06:29.48 


+55:45:45.0 


22.2 x 29.6 


3.54 


29.4 


22.9 


0.41 


0.59 


NGC6946 


Hn 


6.8 


20:34:52.23 


+60:09:14.4 


44.4 X 25.9 


1.25 


120.7 


41.2 


0.11 


39.07 


NGC7331 


LINER 


14.5 


22:37:04.15 


+34:24:55.3 


55.5 x 12.9 


3.56 


-31.1 


26.5 


0.05 


1.22 


NGC7552 


LINER 


21.0 


23:16:10.83 


-42:35:05.5 


37.0 x 29.6 


11.36 


-30.1 


537.4 


0.58 


12.97 


NGC7793 


Hn 


3.8 


23:57:49.84 


-32:35:27.1 


53.6 x 31.4 


0.58 


165.5 


1.0 


0.04 


1.08 



Note. — Col. (2): Based on the optical classification methodology of Kauffmann et al. (2003) applied to the compiled 
spectroscopic sample of Moustakas et al., (2006, in prep.). See also Dale et al. (2006). Col. (3): Adopted distances from 
the literature compilation and local flow model of Masters, et al. (2006, in prep.). Col. (4): Center of extraction in J2000 
coordinates. Col. (6): Physical extraction aperture area. Col. (7): Position angle in degrees E of N of first (slit-aligned) 
axis of extraction rectangle. Col. (8): Total infrared luminosity (3-1100 /im) within the extraction aperture. Col. (9): 
Fraction of global TIR contained in the aperture. Col. (10): Total infrared to B-band luminosity ratio in aperture. 
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the 38" diffraction beam full width of the MIPS 160 /im 
instrument, required a method to overcome the inher- 
ent resolution limits of the two long-wavelength chan- 
nels. We adopted a consistent approach to estimate the 
70 /im and 160 /im surface brightness. Representative 
24/im:70 /tm:160 fim band ratios were obtained in 25" 
radius regions centered on the spectral extraction aper- 
tures, from 24 /im and 70 /mi images convolved to the 
MIPS 160 /im resolution, and from 160 /xm images at na- 
tive resolution. These central color ratios were assumed 
to hold in the smaller spectral apertures, and were com- 
bined with matched 24 /im photometry (which does not 
suffer resolution effects) to predict absolute 70 /im and 
160 /mi surface brightness. From these, the matched to- 
tal infrared intensity, Itir, was calculated according to 
the prescription of Dale & Helou (2002). Except where 
otherwise noted, all references to the total infrared lumi- 
nosity are matched in this way to the spectral apertures. 

Since the 70/im/24/im, and 160/im/24/tm ratios in 
galaxies tend to increase as the strength of the radiation 
field diminishes (e.g. Li & Draine 2001), this method can 
overestimate -Ztir if the innermost regions have warmer 
FIR colors than the surrounding material. To estimate 
the magnitude of this effect, we considered the nuclear 
and inner disk MIPS photometry of Dale ct al. (2005), 
from nearby SINGS galaxies M81 and M51. If the lumi- 
nosity weighted ratio of nuclear to off-nuclear emission 
in the 25" beam is 3:1, the maximum amount by which 
Itir is overestimated using this technique is 25%; most 
galaxies should show no bias in this estimate. 

3.3. IRS Data Reduction and Spectra 

All IRS data were reduced with Cubism (Smith et al., 
2006, in prep.), a custom tool created for the assembly 
and analysis of spatially-resolved spectral cubes from IRS 
spectral maps. The flux calibration was achieved us- 
ing appropriate extended-source corrections (see Smith 
et al., 2006, in prep). All spectral data were processed 
with IRS pipeline version S14. Spectra were extracted 
from IRS modules long- low (LL; 15-37 (im) and short- 
low (SL; 5-15 fim) in matched extraction apertures of 
varying size, chosen to encompass the largest useful re- 
gion of circumnuclear and inner disk emission, including 
star-forming rings and central star formation knots. De- 
tails of the extraction apertures are given in Table 1. 
The enclosed fraction of the total infrared emission il- 
lustrates the range of aperture coverage — from 5-10% 
in large nearby galaxies, to nearly the entire system 
for smaller and more distant targets. Each spectrum is 
composed of data from the four low-resolution orders of 
IRS: SL2 (5.25/im-7.6/jm), SL1 (7.5/im-14.5/im), LL2 
(14.5/mi-20.75/im), and LL1 (20.5 /im-38.5 /im). Slight 
mismatches between order segments, resulting from small 
residual photometric and astrometric uncertainties, were 
addressed by fitting the continuum near the overlap re- 
gion with low-order polynomials, then scaling SL2 to 
match SL1, LL1 to match LL2, and finally SL to match 
LL — a method based on cross-calibration with IRAC 
and MIPS. The scaled spectra were then interpolated 
onto a common wavelength grid and averaged in the over- 
lap region. Typical scaling adjustments were 10%. 

Statistical errors are propagated from the ramp-fitting 
uncertainties produced by the IRS pipeline. Full flux un- 
certainty cubes are extracted alongside the flux cubes, 



and the resulting spectral uncertainties are used con- 
sistently throughout the analysis. At signal-to-noise 
above ~10, systematic errors exceed the statistical un- 
certainties in background-subtracted IRS data, partic- 
ularly in certain wavelength regimes, e.g. 20-22 /mi, 
where LL1 residual fringing is strong, and at the ends of 
the spectral orders, where throughput declines sharply. 
The magnitude of these systematic errors is typically 
15% or less. All spectral mapping BCD exposures were 
background-subtracted, using background frames taken 
from dedicated offset observations or from the periphery 
of the spectral maps themselves. In a small number of 
cases, when suitable SL background observations were 
not available, off-source exposures from the archive, ob- 
served nearby- in-time (At < 3 days) and ecliptic latitude 
(\l — lo \ < 10°), were substituted. 

All 59 low- resolution IRS spectra are shown in Fig. 1. 
The error bars include only statistical uncertainty, 
though fixed uncertainties are important at S/N> 10. 

4. SPECTRAL DECOMPOSITION OF THE MIR 
SPECTRA OF GALAXIES 

In star-forming galaxies, the mid-infrared spectral re- 
gion hosts among the greatest diversity of emission com- 
ponents of any wavelength regime. The broad, strong 
PAH features are set against a background continuum 
transitioning from the rapidly declining direct photo- 
spheric emission of old, cool stars at the shorter mid- 
infrared wavelengths, to intermediate wavelengths dom- 
inated by hot, stochastically-heated small dust grains, 
settling finally on the steadily rising emission of the ener- 
getically important cool dust component, radiating ther- 
mally at longer wavelengths. Set against this diversity 
of continuum and emission features are broad absorption 
bands of silicate (at 9.7 and 18 /im) and, in extreme cases, 
water and other ices (e.g. Spoon et al. 2004). Many of 
the PAH features are blended with prominent forbidden 
spectral lines, with [Neil] 12.8 /mi, blended with a band 
at 12.7/im, being the most severe example. Especially 
at low spectral resolution, separating the lines, features, 
multiple continuum components, and absorption features 
in a consistent way requires a coherent approach. 

To decompose low-resolution IRS spectra, we construct 
a model consisting of a small number of physically moti- 
vated components: dust continuum in fixed equilibrium 
temperature bins, starlight, prominent emission lines, in- 
dividual and blended PAH emission bands, and extinc- 
tion by silicate grains mixed uniformly with all emitting 
components. The main purpose of this simple yet phys- 
ically realistic decomposition method is to estimate ro- 
bustly the strengths of weak features, features blended 
with line emission, and features affected by moderate sil- 
icate absorption — both in absolute terms, and relative 
to the local continuum. 

4.1. The Model 

The surface brightness in the spectral extraction aper- 
ture is modeled as a sum of (1) starlight continuum; (2) 
featureless thermal dust continuum; (3) pure rotational 
lines of H 2 ; (4) fine-structure lines; and (5) dust emis- 
sion features. The model also allows for dust extinction, 
characterized by an optical depth t\ through the emit- 
ting region. If the absorbing dust is assumed to be uni- 
formly mixed with the emitting material, the emergent 
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flux intensity is then: 



Iu = 



t+B v {T+) + 
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£j r (i/) 
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where B v is the blackbody function, is the tempera- 
ture of the stellar continuum, T m are the n = M thermal 
dust continuum temperatures, the I r {v) are the n = R 
resolved (dust) and unresolved (line) emission features, 
and T\ is the dust opacity. These features are described 
in detail below. 



4.1.1. Starlight 

The infrared emission from older stellar populations 
is represented by blackbody emission at fixed tempera- 
ture T* = 5000K, similar to the temperature of the stars 
which dominate stellar emission redward of 3-5 fj,m. This 
temperature was found to best represent the stellar con- 
tinuum in average Starburst99 spectral energy distribu- 
tions (SEDs, Leitherer et al. 1999) for stellar populations 
older than lOOMyr, and wavelength longward of 3 fjja. 
For shorter rest-frame wavelengths, or in fits involving 
near-infrared (NIR) photometry, including a more realis- 
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tic suite of stellar SEDs would be appropriate, as would 
considering the very hot (10 3 K) small grain NIR dust 
emission component inferred in ISO-Phot spectra of nor- 
mal galaxies (Lu et al. 2003). These components are 
unimportant beyond rest frame wavelengths of 3-5 /im. 

4.1.2. Dust Continuum 

We allow up to M = 8 thermal dust continuum compo- 
nents represented by modified blackbodies at fixed tem- 
peratures T m = {35, 40, 50, 65, 90, 135, 200, 300K}, with 
v emissivity normalized at Ao = 9.7 /im. Although 
the bulk of the dust in galaxies is heated to relatively 



low equilibrium temperatures (~15-20K; see Draine et 
al., 2006, in prep.), these grains are too cold to make 
any contribution at A < 40 fim, and are not considered 
here. Contributions from a hot grain continuum, which 
would be important in strong AGN sources, are also not 
included; if present, e.g. for strong SINGS AGN like 
NGC 1566 and NGC5194, these continuum components 
would be absorbed by a combination of warm thermal 
and stellar continuum. The dust continuum is thus de- 
termined by the M parameters r m , each required to be 
non- negative (but not non-zero). The number of com- 



PAH Emission in Star-Forming Galaxies 



7 




30 5 
Rest Wavelength (|a,m) 



10 15 20 30 



Fig. l. 



continued 



ponents was chosen to reproduce smooth distributions 
of grain temperature, though reducing the number (but 
not temperature range) of these components would not 
significantly impact the resulting decomposition. 

No special significance is attached to the strength of in- 
dividual temperature components; rather, they are cho- 
sen to produce a physically realistic dust continuum un- 
derlying the discrete lines and dust features. Includ- 
ing 70 /zm and 160 /um MIPs photometry, along with 
IRAC and near-infrared photometry in matched aper- 
tures, would permit placing more significance on the 



resulting non-stochastic dust temperature distribution; 
this is not pursued here. For extracting PAH, silicate, 
and line strengths, such additional short and long wave- 
length constraints have little impact, and changes in the 
number and exact placement of the temperature bins 
also do not strongly affect the resulting feature param- 
eters. On average, only 5-6 of the maximum allowed 8 
thermal dust components have non-negligible power in 
a given model fit, though all components are well popu- 
lated when considering the sample as a whole. 
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4.1.3. Spectral Line Features 

A diverse collection of spectral lines is emitted in the 
mid- infrared range probed by the IRS, including the 
H2 S(0)-H2 S(6) pure rotational lines of molecular hy- 
drogen, strong H 11-region and photo-dissociation region 
(PDR) cooling lines from low-ionization species such as 
Si + , with an ionization potential of 8.2 eV, to massive 
star and AGN tracers from 0+++ (55eV) and Ne 4+ (97 
cV). These lines are modeled primarily to allow realis- 
tic decomposition of blended features and the underlying 
continuum. The strong line pair ratio [Nem]/[Neii] pro- 



vides an indicator of radiation hardness unaffected by all 
but the heaviest dust extinction (Thornley et al. 2000). 

Spectral lines in the model are represented by Gaussian 
profiles, with widths set per instrument module based on 
direct measurements of unresolved features, and shifted 
into the rest frame. The lines which are fitted to the spec- 
tra are listed in Table 2. Line wavelengths are allowed to 
vary by 0.05 fim (roughly 1/2 of a spectral pixel, approxi- 
mately the wavelength calibration uncertainty at the blue 
ends of each IRS order), and line widths are allowed to 
vary by 10% from their initial default values. The central 
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intensity of all lines is required to be non-negative. 

4.1.4. Dust Features 

Dust features are represented by individual and 
blended Drudc profiles: 



j(r) 



(A/A r -A r /A) 2 + 7r 2 



(2) 



where A r is the central wavelength of the feature, j r is the 
fractional FWHM, and b r is the central intensity, which is 
required to be non-negative. The Drude profile is the the- 



oretical frequency profile for a classical damped harmonic 
oscillator, and is a natural choice to model PAH emission 
features (and other emission and absorption bands, e.g., 
the 2175A extinction feature, Fitzpatrick & Massa 1986). 
The Drude profile, like the Lorentzian, has more power 
in the extended wings than a Gaussian. Power in the 
overlapping wings of multiple bands can dominate the 
apparent continuum, especially in the 5-12 /tm region. 
The integrated intensity of the Drude profile is 



j(r) 



r ] dv 



7TC b r J r 
IT 



(3) 



10 



Smith et al. 



TABLE 2 
Unresolved Line Parameters 



Line 


A r ( /im) 


FWHM(ttm) 


(1) 


(2) 


(3) 


H2 S(7) 


5.511 


0.053 


H2 S(6) 


6.109 


0.053 


H 2 S(5) 


6.909 


0.053 


[Aril] 


6.985 


0.053 


H 2 S(4) 


8.026 


0.100 


[Ar III] 


8.991 


0.100 


£12 0(0) 


y.Doo 


U. 1UU 


[Siv] 


10.511 


0.100 


H 2 S(2) 


12.278 


0.100 


[Nell] 


12.813 


0.100 


[Nem] 


15.555 


0.140 


H 2 S(1) 


17.035 


0.140 


[Sin] 


18.713 


0.140 


[Oiv] 


25.910 


0.340 


[Fen] 


25.989 


0.340 


H 2 S(0) 


28.221 


0.340 


[Sin] 


33.480 


0.340 


[Sin] 


34.815 


0.340 


Note. 


- Col. (1 


) : line name. 



Col. (2): central wavelength. 
Col. (3) : full- width at half max- 
imum. 



Table 3 lists all modeled dust features, some of which 
are represented by up to 3 sub-features, grouped as 
blended complexes. The 17 /im complex has two re- 
solvable components, at 16.4 /im (Moutou et al. 2000) 
and 17.4 /im (Smith et al. 2004; Werner et al. 2004) - 
see Fig. 4. For the other blends (7.7/im, 11.3 /im, and 
12.7 /im), sub-features were adopted to best reproduce 
the band shape. Unresolved complexes of this type are 
considered only in aggregate, and the sub-features com- 
prising them are not intended to represent distinct phys- 
ical emission bands. The number and placement of these 
sub-features was determined by training the model on a 
subset of 25 high signal galaxies from the sample, by fix- 
ing all other fitting components, and allowing the width 
and central wavelengths to vary. A similar optimization 
was performed for the weaker individual features. Once 
determined, feature positions and widths were held fixed. 
Related fitting approaches have been used with ISO spec- 
tra (Boulanger et al. 1998; Verstraete et al. 2001). 

Features from 5-10 [im — The bulk of the PAH emis- 
sion between 5-12 /im is comprised of the well-studied 
set of features discovered in planetary nebulae by Gillctt 
et al. (1973), and probed in detail in Galactic objects 
with ISO (e.g. reflection nebula NGC 7023, Cesarsky 
et al. 1996). The weak 8.33 /im band, which sits between 
the prominent 7.7/im and 8.6 /im features, was tenta- 
tively detected in the ISO SWS spectra of M82 (Sturm 
et al. 2000) , and is present in many of our spectra, serv- 
ing in most cases only to fill in the saddle between the 
two stronger features. A weak new feature in external 
galaxies occurs at 6.69 /im (see also Werner et al. 2004). 

A small number of stellar-dominated spectra from el- 
liptical and other early type systems in the sample show 
broad 10 /im excess (e.g. NGC 4552), which has been 
attributed in similar systems to silicate emission from 
the windos of AGB populations (Bressan et al. 2006). 
These features are not considered in our fit, and are weak 



TABLE 3 

Dust Feature Parameters 



A r ( (im) 


7r 


FWHMOm) 


(1) 


(2) 


(3) 


^ 97 


U.U04: 


n 1 7Q 


^ 7ft 


u.uoo 


u.zuu 


6 22 


u.uou 


n 1 &7 

u. 10 1 


u.uy 


n n7n 




7 A9 a 


n 1 o(\ 
u. 1ZO 


u.yoo 


/ .OU 


n C\AA 




7 

* .00 


n 

v.yjoo 




8.33 


0.050 


0.417 


8.61 


0.039 


0.336 


10.68 


0.020 


0.214 


11.23 b 


0.012 


0.135 


11.33 b 


0.032 


0.363 


11.99 


0.045 


0.540 


12.62 c 


0.042 


0.530 


12.69 c 


0.013 


0.165 


lo.4o 


U.U4U 


U.ooy 


14.04 


0.016 


0.225 


14.19 


0.025 


0.355 


15.90 


0.020 


0.318 


16.45 d 


0.014 


0.230 


17.04 d 


0.065 


1.108 


17.375 d 


0.012 


0.209 


17.87 d 


0.016 


0.286 


18.92 


0.019 


0.359 


33.10 


0.050 


1.655 


Note. 


- Col. 


(1): central 



wavelength. Col. (2): frac- 
tional FWHM. Col. (3): full- 
width at half maximum. 

a 7.7/im Complex 

b 11.3/im Complex 

c 12.7/im Complex 

d 17/tm Complex 



enough not to influence the recovery of stellar continuum 
in these cases, though they can induce spurious dust fea- 
tures (especially at 7.7/im), and for this reason are ex- 
cluded from results pertaining to these features. 

Features from 10-20 /im — Along with the strong bands 
at 11.3 /im, 12.7/im, and 17.0/im, a variety of weak fea- 
tures are seen in most low-resolution spectra between 
10-20 /im. The 10.7/im band, which falls in the (po- 
tentially silicate-extincted) trough between the main 7.7 
and 11.3 /an features was tentatively identified in the ISO 
spectra of M82 (Sturm et al. 2000), and is not uncom- 
mon in the sample, occasionally blending with nearby 
[Siv] at 10.51 /mi. Other weak bands arc found at 
12.0 /im, 12.6 /xm, 13.6 /mi, 14.3 /xm, 15.9 /im, 16.4 /im, 
and 17.4 /im. Sturm et al. (2000) identified all of these 
features in ISO SWS spectra of the bright starburst 
galaxies NGC 253 and M82, and many of the individ- 
ual sub-features were first discovered in Galactic sources 
(e.g. Moutou et al. 2000). The 14.3 /im dust feature may 
be blended with [Nev] 14.32 /im in sources with AGN 
nuclei. Another very weak feature at 18.92 /im is likely 
the same as the 19 /im feature found in Spitzer spectra of 
reflection nebula NGC 7023 by Werner et al. (2004). It 
is blended with [S ill] at 18.7/im. Dust features >10/im 
are marked in Fig. 2 on the high-resolution IRS nuclear 
spectrum of NGC 5195 (M51b). Note in particular that 
the 12.7/im PAH band, which is blended with [Neil] at 
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Fig. 2. — High-resolution spectrum of NGC5195 (M51b) illustrating the sub-set of features from 10— 19 /tm identified in the full low- 
resolution sample. The PAH feature central wavelength and FWHM for features from Table 3 with central wavelengths >10/^m are 
indicated above the spectra, and unresolved emission lines in the spectrum are also labeled. 



low-resolution, cleanly separates with [Neil] to the red 
wing of the blend at higher resolution. 

Interestingly, the strongest feature at A > 12 /im, the 
broad band centered at 17.04 /im, was not seen in the 
ISO spectra of the brightest starbursts, and was first 
identified in an extragalactic source in the SINGS galaxy 
NGC7331 (Smith et al. 2004). Despite the strength of 
the 17 /mi complex, which can carry power comparable 
to the prominent 11.3 /mi band (see § 5.4), its large width 
(~1.1 /im), and the lack of available low- resolution spec- 
troscopy at wavelengths above 16 /im combined to make 
this prominent band difficult to detect with ISO. The 
17 /mi feature was reported simultaneously in Spitzer 
spectra of reflection nebulae NGC 7023 (Werner et al. 
2004). Evidence for a similar band complex was found 
in ISO SWS spectra of Galactic sources (Beintema et al. 
1996; Van Kerckhoven et al. 2000). 

Features beyond 20 ' [im — A single feature was in- 
cluded beyond 20 /mi: a broad band centered at 33.1 /mi. 
Though not present in all sources, it underlies the [Si n] 
and [S in] lines, and is required to obtain accurate line 
fluxes for these lines. Sturm et al. (2000) notice a broad 
plateau from 33-34 /im in 3 ISO SWS galaxy spectra, and 
attribute it to crystalline silicates. A similar plateau was 
noted in the SWS spectra of planetary nebulae (Waters 
et al. 1998). This feature is strong in ~l/6 of our sample 
(e.g., see NGC 4321), with no correlation of the feature's 
strength with silicate extinction (see §5.2). 

4.1.5. Dust Extinction 

We allow for absorption by dust with extinction prop- 
erties similar to Miky Way dust. We take the infrared 
extinction to consist of a power-law plus silicate features 
peaking at 9.7 /im and 18 /mi: 



r(A) =r 9 . 7 [(l-/3)Psi(A) + /3(9.7/ i m/A) 1 - 7 ] (4) 
where T9.7 is the total extinction at 9.7 /im, and 
Psi(A) is the normalized silicate profile function, with 
Psi(9.7/mi) = 1. We adopt (3 — 0.1; with this choice, 
r(3.6/im) = 0.54 r(9.7 /mi), which is within the cur- 
rent uncertainty range regarding extinction from 3-8 /tm 
(see, e.g., Draine & Li 2001). The 9.7/mi component 
of -Psi(A) is taken from the Galactic center profile of 
Kemper et al. (2004) between 8.0 /im and 12.5 /tm, with 
Psi(8/im) = 0.06. For A < 8 /im, we smoothly extend 
this profile by assuming P Si (A) = P si (8 /im) e 2 03 ( A " 8 ^ m ) . 
For A > 12.7 /tm, we add a Drude profile extinction 
component peaking at 18 /im, with FWHM=4.45 /tm 
(7 r =0.247), smoothly joining the two profiles between 
12.5 /tm and 12.7/im. The fitted silicate depth T9.7 is 
required to be non-negative; no silicate emission was in- 
cluded in the model. We also explored using the Galactic 
center extinction curve of Chiar & Tielens (2006), with 
similar values of T9.7 derived, but the much larger ex- 
tinction at long wavelengths {A^Q^ m ~ As^m) produced 
poorer overall fits in absorbed glaxies with our choice of 
dust continuum temperatures. The adopted extinction 
curve is shown in Fig. 3. 

The default extinction geometry is taken to be absorb- 
ing dust fully mixed with emitting stars and grains; a uni- 
form foreground screen of dust is also considered. Since 
all components of the model are presumed to be equally 
affected by dust extinction, line and feature strengths are 
automatically corrected to their extinction- free values. 
This implicit correction may be unphysical if the material 
emitting the mid-infared continuum, from which the dust 
absorption is primarily determined, and the PAH emit- 
ting material undergo differential extinction. The TIR 
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Fig. 3. — The adopted extinction curve, along with the curve of 
Chiar & Tielens (2006). 

luminosity estimates described in § 3 are not extinction- 
corrected, but the correction for A > 25 /im, where the 
bulk of the TIR power emerges, is negligible. 

4.2. Fitting and Uncertainties 

The model in Eq. 1 is fitted to the observed flux inten- 
sity by minimizing the global % 2 using the Levenberg- 
Marquardt algorithm, modified to allow the values of 
individual fit parameters to be fixed within allowed lim- 
its. Formal statistical uncertainties, derived from the IRS 
pipeline ramp fitting uncertainties, are returned for each 
fit parameter, and full uncertainties in combined quan- 
tities such as integrated intensities are formulated using 
the full covariance matrix. In particular, for estimat- 
ing uncertainties in the integrated intensity of blended 
dust features, correlations between components can in- 
troduce significant cross-term corrections to the propa- 
gation of error. For upper limits, additional uncertainty 
estimates, based on the residuals of the observed flux in- 
tensity from the best fitting model, are adopted. The 
response of the model to a variety of simulated changes 
to the continuum was tested, and the strengths of dust 
features recovered were found to be robust against such 
variations. In particular, adding strong stellar contin- 
uum to PAH-dominated spectra and re-fitting resulted 
in feature fluxes which changed by less than 2%. 

4.3. Example Decomposition and Tool Availability 

The detailed components of the spectral decomposi- 
tion method are illustrated in Fig. 4 for NGC 1482, and 
in Fig. 5 for four additional examples — two sources 
with significant dust extinction, and two sources with 
starlight dominating their short-wavelength continuum. 
All galaxies in the sample are well reproduced by the 
modeled components. A tool implementing the decom- 
position method is available in Pahfit, an IDL package 
developed for the decomposition of low resolution IRS 
spectra comprised of stellar and thermal dust continuum 
(not including very hot dust in strong AGN), dust fea- 
tures, line emission, and silicate absorption. Source code 
and documentation for Pahfit are available 14 . 

5. RESULTS 

The spectral decomposition of all 59 sample galaxies 
was used to investigate silicate extinction, the energetics 
of PAH emission, variations in the large number of indi- 
vidual and blended PAH features recovered, and the im- 
pact of these variations on studies of high-redshift galax- 
ies in deep broadband surveys. Recovered strengths for 

14 http://sings.stsci.edu/pahfit 
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all important dust features, and the [Neil] and [Nem] 
lines are given in Table 4. 

5.1. The Nature of the 11 [im Emission Band 

The 17/im PAH complex is the strongest band long- 
ward of 12 ^m. In contrast to the highly variable plateau 
reported by Peeters et al. (2004), all SINGS spectra ex- 
hibit a dominant component of the 17 /im complex which 
varies very little in position and width. This unifor- 
mity could be a result of beam-averaging over a range 
of physical conditions in PDRs drawn from similar en- 
vironmental distributions, or perhaps an indication that 
the Galactic sources considered by Peeters et al. were 
dominated by atypical band carriers or radiation envi- 
ronments. The distinct flanking sub-features at 16.4 /im 
and 17.4 /im contribute only ~20% of the total power in 
the complex, with the remainder contained in the broad 
(FWHM'vl /jm) feature centered at 17.04 /im. A small 
number of galaxy spectra in the sample do exhibit an 
apparent plateau-like extension of the 17/im band to 
~19/im (e.g. NGC 1566, NGC 4569), but it does not 
correlate with the strength of the 17/im feature or any 
other PAH band, and may arise from a different carrier. 

Although laboratory spectra of PAH mixtures contain 
many features between 15-20 /im presumed due to C- 
C-C out-of-plane bending modes (Peeters et al. 2004), 
it is of interest to compare the general behavior of the 
17/im band with the other more traditional PAH bands. 
Fig. 6 compares the integrated feature power of the 17/im 
and 11.3 /im complexes. The features correlate very well 
with each other over 2 orders of magnitude, strongly sup- 
porting the association of the 17 /im complex with PAH 
emission. Some variations in the 17/im/ll /im intensities 
ratio is evident; these will be explored in § 5.6. 

5.2. The Incidence of Silicate Extinction 

Helou et al. (2000) found little evidence for silicate 
absorption in the 5.9-11.7/im spectra of the ISO Key 
Project sample of 45 normal galaxies, and Sturm et al. 
(2000) also concluded that high-resolution ISO spectra 
of bright starbursts like M82 were consistent with little 
or no silicates. With a more varied sample, and access to 
the full 9.7 /im and associated 18 /im features, the SINGS 
spectra are uniquely suited to the investigation of mod- 
erate silicate absorption in a range of galaxy types. 

Attempting to measure weak to moderate silicate ab- 
sorption in spectra with strong PAH emission features 
can lead to ambiguous measures of both the feature 
strengths and the silicate optical depth. In systems dom- 
inated by deeply embedded infrared sources (e.g. many 
ULIRGs) with strong silicate absorption, strong to weak 
or even absent PAH bands, and accompanying shorter 
wavelength ice absorption features, the identification and 
measurement of the silicate and other absorption optical 
depths is straightforward (e.g. Armus et al. 2004; Spoon 
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Fig. 4. — Detailed decomposition of NGC 1482 from 5—37 /im. Red solid lines represent the thermal dust continuum components, the 
magenta line shows the stellar continuum (weak in this case), and the thick gray line shows the total (dust + stellar) continuum. Blue 
lines set above the total continuum are dust features, while the narrower violet peaks are unresolved atomic and molecular spectral lines. 
All components are diminished by the fully mixed extinction, indicated by the dotted black line, with axis at right. The solid green line is 
the full fitted model, plotted on the observed flux intensities and uncertainties. 



et al. 2005; Roussel et al. 2006). However, as the sili- 
cate optical depth decreases, and the PAH strength rela- 
tive to the local continuum increases, the PAH bands at 
7.7/8.6/im and 11.3 /im, which closely flank the main sil- 
icate absorption trough centered at 9.7/xm, can mask or 
mimic the true silicate absorption. With this absorption 
feature alone it is difficult to discriminate between mod- 
erate PAH emission superposed on a silicate-absorbed 
continuum, and strong PAH features with a relatively 
weak underlying continuum. However, since the dust 
continuum at 15-20 /im is fairly regular in the spectra 
of star-forming galaxies, increasing silicate absorption in 
the associated 18 /im band results in a measurable change 
in the shape of the 15-20 /im continuum. 

A compelling example of this degeneracy in recover- 
ing silicate extinction from limited wavelengths around 
the 10 /im band is found NGC 5866. Lu et al. (2003) 
noted an unusually high value of L(6.2/tm)/L(7.7/im) 
for this compact, edge-on disk galaxy. With access to 
the full 5-38 /im spectrum, and identification of the non- 
negligible silicate extinction 79.7 = 0.36, the band ratio 
L(6.2 /xm)/L(7.7 /im) is placed on the low end of the dis- 
tribution (see Fig. 13). In this case, the strong pho- 
tospheric continuum which dominates NGC 5866 's spec- 
trum out to 10 /im may have further complicated accu- 
rate measurement of the band strengths. 

Only 13.6%, or 8 of the 59 sample spectra shown in 
Fig. 1, require non-negligible absorption by silicate dust 



TABLE 5 
Silicate Dust Extinction Sources 





T9.7 






Galaxy 


Mixed 


Screen 


X 2 (r S i = 0)/x 2 (T Si ) 


NGC 1266 


2.05 


0.91 


1.868 


NGC 1482 


1.04 


0.37 


1.346 


NGC 3198 


4.88 


1.90 


4.978 


NGC 4536 


0.37 


0.18 


1.086 


NGC 4631 


0.40 


0.20 


1.031 


NGC 5033 


0.41 


0.20 


1.058 


NGC 5866 


0.36 


0.17 


1.060 


NGC 6946 


0.76 


0.37 


1.258 



to reproduce their spectra. These objects, and their in- 
ferred silicate optical depths 79.7, are listed in Table 5, for 
the default fully mixed dust geometry, as well as a uni- 
form foreground screen. The silicate depths are lower 
for the screen geometry, by roughly a factor of two. 
The mixed and screen models provide equally accept- 
able fits. Also given is X 2 ( r Si = 0)/x 2 (Tsi) — the ratio of 
X 2 ( T Si = 0), formed by explicitly forcing the silicate depth 
to be zero, to the normal X 2 ( T Si) including absorption. 
This ratio demonstrates the degree to which the fit is im- 
proved by adding silicate extinction, and indicates that 
the decomposition method can robustly recover modest 
extinction with mixed silicate optical depths rsi > 0.35, 
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Rest Wavelength (urn) 

Fig. 5. — Example decomposition of four SINGS spectra, for sources with and without significant silicate extinction, and dominated by 
either stellar or thermal dust continuum at short wavelengths. All plot symbols as in Fig. 4. Logarithmic wavelength axis. 
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Fig. 6. — The correlation between the intensity of the 11.3 /im 
and 17 (im PAH complexes, shown with a slope unity line. 



or screen depths half of that value. Thus it appears 
that most normal galaxies do not exhibit dust extinc- 
tion to the mid-infrared emitting material in their cen- 
tral few square kiloparsecs above Ay ~ 5 for fully mixed 
dust geometries, or Ay ~ 3 for foreground screens of 
dust (assuming A v /t S [ — 16.6, Rieke & Lebofsky 1985). 
Screen extinctions derived separately from optical nu- 
clear spectra for the SINGS sample show (Ay) ~ 1.3, 
consistent with the inferred silicate depths (Dale et al. 
2006). Brandl et al. (2006) found a larger range and 
greater incidence of silicate absorption depths in the nu- 
clear spectra of a sample of starburst galaxies. 



5.3. Recovering PAH Strengths 

A widely used method for extracting PAH strength 
involves selecting narrow continuum regions free from 
line or feature emission, forming linear or cubic spline 
fits through these continuum points, and integrating the 
spectrum above this pseudo-continuum (e.g. Uchida et al. 
2000; Peeters et al. 2002; Brandl et al. 2006). An applica- 
tion of this method to the SINGS spectrum of NGC 2798 
is shown in Fig. 7. The method's advantage is that it is 
relatively straightforward to apply, and requires no spe- 
cific assumptions regarding the nature of the continuum. 
The disadvantage is that locating nearby continuum an- 
chors free from line, feature, or absorption contamination 
can prove difficult and to some degree arbitrary, with the 
resulting feature strengths depending sensitively on the 
exact choice of fiducial continuum points. This problem 
is compounded in low quality spectra, where continuum 
regions free from these contaminating effects may be less 
readily identified. 

To compare the band strengths recovered by the spline- 
fitting and full decomposition techniques, we have ap- 
plied the cubic spline method for the features at 6.2 /xm, 
7.7/xm, 8.6 /im, and 11.3 /im to the full sample, fixing 
the rest wavelength continuum pivots following Peeters 
et al. (2002). Figure 7 illustrates the chosen fixed rest- 
frame continuum anchor points: 3 points for 6.2 /zm, 4 
for 11.3 /im, and 5 for the 7.7+8.6 /mr blend. The PAH 
bands at 12.6 /im and 17/zm are blended with strong 
emission lines ([Neil] and H2S(1)), and for this reason 
cannot be recovered with spline-based methods. 
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Wavelength (|im) 

Fig. 7. — An example comparing the spline-based and full feature+continuum decomposition in SINGS galaxy NGC 2798. Three different 
spline fits are shown with continuum anchor points: in green (leftmost), underlying the the 6.2 /an feature, in red (center), underlying both 
the 7.7 /an and 8.6 /an features, and in orange (rightmost), under the 11.3 /im feature. Vertical dotted lines mark the ranges over which 
the spline-subtracted features are integrated. The fitted continuum, dust, and line features, are as in Fig. 4. 




TABLE 6 

PAH Feature Luminosity Estimates: 
Full Decomposition vs. Spline Ratios 



10.00.1 1.0 
PAH Equivalent Width ((im) 

Fig. 8. — The ratio between the PAH features strength from the 
full spectral decomposition, to the PAH strength recovered from 
spline-based continuum fitting vs. the equivalent width measured 
from the decomposition. Mean and unity values are indicated with 
horizontal dotted lines. 



Fig. 8 shows the PAH strengths recovered from full 
feature decomposition vs. spline fitting as a func- 
tion of the feature equivalent width (defined as J [I v — 
I^ ont ) / Iy° nt dX). Equivalent widths of the shortest wave- 
length features can be undefined, since the underly- 
ing continuum can vanish. In these cases, the profile- 
weighted average continuum was substituted. For the 
two unblended features in clean continuum regions, at 
6.2 /im and 11.3 /im, the decomposition and spline-based 
methods yield feature powers differing by an approxi- 
mately constant ratio, L(decomp)/i (spline) ~1.75, in- 
dependent of the equivalent width of the feature. Since 
a fraction f W i ng ~ 59% of the power in a Drude profile 
lies outside its FWHM, it is not surprising that the offset 
between the two is roughly 1/ f W i ng . 

For the blended features at 7.7 /im and 8.6 /im, how- 



Feature 


Z/(decomp)/L(spline) 


6.2 fim 


1.70 ±0.42 


7.7 /im Complex 


3.53 ±0.62 


8.6 /tm 


3.31 ±0.64 


11.3 fim Complex 


1.86 ±0.15 



ever, the fraction of the full feature strength recovered 
by the spline method is lower, with larger scatter. This 
results from the continuum anchor placement in the 
filled trough between these features, located to ensure 
that the de-blended peaks are independent. The over- 
all scatter is probably somewhat larger than would have 
been achieved with hand-tuned spline continuum place- 
ment in all spectra. However, the large systematic off- 
set would remain, such that the spline-based methods 
underestimate by a factor of 3-6 the full 7.7 /im and 
8.6 /im feature strengths. Table 6 lists the average off- 
sets for feature equivalent widths larger than 0.1 /im, af- 
ter applying an iterative 3tr trimming of outliers. These 
L(decomp)/L (spline) factors can be used to scale spline- 
recovered PAH fluxes to estimate crudely the total power 
in the PAH features - although, if possible, a more ac- 
curate estimate will be obtained directly using a decom- 
position method like the one described in § 4. They may 
have dependence on the form of underlying continuum 
and presence of silicate or other absorption, and thus 
may not be appropriate for all spectra. 

5.4. PAH Energetics 
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Fig. 9. — The fractional PAH luminosity relative to the total PAH power for the main features identified in the spectra, against the 
integrated PAH contribution to the total infrared luminosity from 3-1100 fim, in the same region. Blended complexes are marked by a 
double arrow in the legend; see Table 3. 
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Fig. 10. — The distribution of the integrated PAH luminosity Fig. 11. — The PAH to MIPS 24 /mi luminosity ratio vs. far- 
relative to the total infrared, ipAH/^TIRi f° r the entire sample, infrared color for the four strongest PAH bands, in matched aper- 
and individually for galaxies with AGN and Hn/Starburst classi- tures. 
fications 



The total power emitted in the PAH bands can be 
large, and the current sample provides an opportunity 
to form a census of the energetic importance of individ- 
ual features, independent of the underlying mid-infrared 
continuum, which can be dominated by widely differing 
emission processes and environments and depend on the 
details of the stellar and dust grain populations. In Fig. 
9, the fractional power relative to the total PAH lumi- 
nosity is shown for all significant features present in the 
sample. 

The 7.7 /im PAH complex can contribute nearly one- 
half of the total PAH luminosity, and deliver up to 10% 
of the total infrared luminosity alone. In aggregate, the 
combined power output of all PAH features ranges from 



a few percent up to 20% of the total infrared. Fig. 10 
shows the distribution of ipAH/-^Tm in the apertures 
considered, for the full sample, and separately for galax- 
ies with evidence of nuclear Seyfert or LINER emission 
in their nuclear optical spectra, and galaxies with only 
Hll-region or starburst classifications. The peak is near 
-£tah/-£'TIR = 0.1 for both H II and AGN galaxies. At the 
highest relative PAH luminosities, systems dominated by 
star formation predominate (at £pah/£tir > 0.12, only 
3/14 systems host AGN), despite the fact that the spec- 
tral apertures include both nuclear and inner disk emis- 
sion. This may suggest a natural limit on the absolute 
PAH strength in sources with weak AGN, either due to 
partial destruction of the grains, or the contribution of 
other continuum sources to the total infrared luminosity. 
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Fig. 12. — The distribution among the sample of the fractional 
contribution of dust features to the Spitzcr broad bands, for IRAC 
8 fim at redshift z = 0, with and without stellar continuum removed 
from the band, as well as MIPS 24 (im at z = 2. 



Fig. 11 demonstrates the large scatter in the ratio 
L(PAH)/j/L 1/ (24 /tm) for the four strongest PAH features 
as a function of the local far-infrared color temperature 
vL v (2A /im)/VL„(70 /an). The full range is typically over 
a factor of 10 for each band, and may in fact be a lower 
limit on the true scatter in the global emission of galax- 
ies, for which diffuse scattered light in the outer regions 
can account for a substantial fraction of the 24/tm lu- 
minosity (Calzetti et al. 2005). For all bands, there is 
a definite negative trend of L(PAH)/i/L„(24 fim), which 
suggests that a mix of both small stochastically-heated 
and thermally emitting grains contribute to the 24 /mi 
continuum, perhaps due to varying contributions from 
PDR and Hll regions (e.g. Laurent et al. 2000). 

5.5. The PAH Contribution to IRAC 8 /im and MIPS 

24 [am 

With the full spectral decomposition of continuum and 
dust features, we can examine the fractional contribu- 
tion of PAH emission and stellar or thermal dust con- 
tinuum to the broad Spitzer bandpasses. We consider 
IRAC 8 /mi at z = 0, and MIPS 24/tm at z = 2. These 
broad bands are commonly used as surrogates for PAH 
strength in nearby and distant galaxy surveys (e.g. Hogg 
et al. 2005; Papovich et al. 2006; Reddy et al. 2006). 
Since IRAC 8 /mi at low redshift, and MIPS 24 /mi at 
z <~ 2 encompass the 7.7 /im and 8.6 /im PAH bands — 
typically the strongest (although, see § 5.6) — a common 
assumption is that broadband photometry at these rest 
wavelengths serves as a direct proxy for PAH emission, 
total infrared or bolometric luminosity, and, potentially, 
star formation rate. For low redshift sources, prescrip- 
tions exist for removing the stellar contribution using the 
3-5 /jm channels of IRAC (e.g. Helou et al. 2004). The 
relatively isolated 24 /jm band does not typically permit 
such starlight removal (although the IRS Peak-Up Blue 
channel at 16 /im could in principle serve this purpose). 

Fig. 12 highlights the large fractional contribution of 
PAH feature emission to the uncorrected IRAC 8 /tm 
band at z = and MIPS 24 /im band at z = 2. The 
median contribution for both ^0.7, with a small number 
of sample members (sources dominated by stellar con- 
tinuum) showing very little PAH contribution. For over 
75% of the sample, PAH emission contributes more than 



one half of the broad band 8 /im power. If the stellar con- 
tinuum is perfectly removed, the distribution shifts only 
slightly to higher fractional PAH contribution. 

The typical source in deep MIPS 24 /im surveys has to- 
tal infrared luminosity Ltir > 10 11 L Q (Papovich et al. 
2006; Yan et al. 2004), and may therefore exibit a much 
different contribution of PAH emission in the 8 /tm rest 
frame compared to star-forming galaxies of lower lumi- 
nosity. The sample explored here has median global Ltir 
of 0.8 x 1O 1O L . It should also be noted that the domi- 
nance of PAH emission in the band, e.g. MIPS 24 /tm at 
z^2, does not imply anything about the efficacy of using 
this single photometric point to derive the total infrared 
luminosity or star formation rate; see Dale et al. (2005). 

5.6. PAH Strength Variations 

Fig. 13 demonstrates the variation in inter-band 
strength ratios over the full sample for the four strongest 
PAH bands at 6.2 /tm, 7.7/tm, 11.3 /tm, and 17 /tm. Both 
the 6.2/zm/7.7/im and 17/im/11.3/im ratios vary by 
more than a factor of two. Sources with weak AGN 
are almost entirely offset from sources with Hn-like 
nuclear optical spectra — an effect which will be ex- 
plored in the next section. Sources with non-negligible 
mixed silicate absorption optical depths (circled) are off- 
set slightly from the main locus of points for their type 
(Hll or AGN), to lower values of L(6.2/mi)/L(7.7/mi) 
and L(11.3 /tm)/L(17/im). This is in the direction of in- 
creased reddening, despite the implicit de-reddening per- 
formed by the spectral decomposition method, and could 
be evidence of different absolute extinctions for the un- 
derlying continuum and the PAH emitting regions, or of 
limitations in the adopted extinction curve. 

The medians and 10%-90% range of variation for lu- 
minosity ratios among the most important PAH bands, 
the integrated PAH luminosity, and the total infrared 
luminosity are given in Table 7. Typical variations in 
the inter-band luminosities are a factor of 2-5, with la 
variations in these ratios of ±30%. 

5.6.1. PAHs in AGN 

Nuclear black holes are now recognized as a fundamen- 
tal property of many local elliptical and spiral galax- 
ies: >l/3 have optical nuclear spectra typical of low- 
luminosity active galactic nuclei (Ho et al. 1995, 1997). 
The SINGS sample was not defined to include power- 
ful AGN sources, yet roughly half of the galaxies con- 
sidered here show some form of AGN activity in their 
nuclear optical spectroscopy. The SINGS AGN sources 
are in the low- luminosity regime, with L([0 ill] 5007A) < 
1O 6 5 L (compared to L([Om]) ~ 10 8 -10 9 L Q for tradi- 
tional AGN, Kewley et al., 2006, in prep.). This sample 
makes it possible to test the behavior of PAH emission 
in the presence of weak AGN. 

In Fig. 14, the ratio of two strong PAH bands is 
compared to the hardness indicator [Ne m]/[Ne 11], with 
an ionization potential ratio of 41eV/21.6eV. Galax- 
ies with H 11 region or starburst-likc optical spectra 
exhibit a nearly constant L(7.7 /im)/L(11.3 /im) ratio 
across the full range of radiation hardness, includ- 
ing those low metallicity systems (when PAHs are de- 
tected) with [Ne in] /[Neil] > 2. In striking contrast, 
the AGN as a group are offset below the locus of H II- 
type galaxies, even at moderate hardness ratios, and 
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Fig. 13. — The inter-band strength ratios of the four main PAH bands at 6.2 fim, 7.7 /im, 11.3 (im, and 17 /im. Triangles (red) indicate 
galaxies with AGN (Seyfert or LINER) types, whereas squares (blue) are galaxies with H Il-dominated nuclei. Galaxies with measurable 
silicate absorption (rgj > .01) are circled. Also shown are the extinction trends for these band ratios, parametrized by the silicate optical 
depth for both uniform foreground screen and fully mixed dust geometries. The highlighted regions indicate zones used for the construction 
of template spectra (see §6). 



TABLE 7 

PAH Band Luminosity Ratios L(Ai)/L(A2) 



Band^\( 2 ) 


6.2 pm 


7.7pm** 


8.6 pm 


11.3 |im s 


12.6 |im tt 


17 (im* 


EPAH 


TIR/100 



6.2 pm •■■ 0.28 ( 0.21-0.75) 1.5 ( 1.2-3.0) 1.1 (0.52-1.5 ) 1.8( 1.1-2.8) 1.9(0.93-2.9) 0.11 (0.086-0.15 ) 1.1 (0.35-2.1 ) 

7.7pm 3.6 ( 1.3-4.8 ) ■■■ 5.7 ( 4.7-9.0) 3.6 ( 1.5-4.8 ) 6.5 ( 3.6-8.3) 6.9 ( 3.1-11 ) 0.42 ( 0.18-0.45 ) 4.1 (0.54-6.2 ) 

8.6 pm 0.66(0.38-0.86) 0.18( 0.11-0.21) ■•■ 0.68(0.35-0.81) 1.1 (0.66-1.5) 1.2(0.65-2.0) 0.073 (0.049-0.088) 0.72 (0.16-1.2 ) 

11.3pm 0.95(0.65-1.9 ) 0.28 ( 0.21-0.67) 1.5( 1.2-2.9) ■■■ 1.8 ( 1.5-3.2) 1.9 ( 1.4-2.8) 0.12 (0.092-0.18 ) 1.1 (0.57-1.7 ) 

12.6 pm 0.55 (0.36-0.93) 0.16 ( 0.12-0.27) 0.93(0.65-1.5) 0.57(0.31-0.65) ■■■ 1.1 (0.54-1.6) 0.065(0.047-0.087) 0.60(0.19-0.92) 

17 pm 0.52 (0.34-1.1 ) 0.15 (0.089-0.32) 0.81 (0.50-1.5) 0.53 (0.35-0.72) 0.93 (0.62-1.8) ■■■ 0.059(0.039-0.097) 0.62(0.28-0.94) 

X)PAH 8.8 ( 6.6-12 ) 2.4 ( 2.2-5.7 ) 14 ( 11-20 ) 8.6 ( 5.5-11 ) 15 ( 11-21 ) 17 ( 10.-25 ) ■■■ 10. ( 3.2-16 ) 

TIR/100 0.88(0.48-2.9 ) 0.25 ( 0.16-1.9 ) 1.4(0.85-6.1) 0.88 (0.60-1.7 ) 1.7 ( 1.1-5.2)1. 6 ( 1.1-3.6) 0.100(0.064-0.31 ) ■__ 

Note. — Each entry is encoded as med (low-high) where med is the median of the ratio, and low-high is the 10%-90% range of variation. 

''Indicates blended PAH Complexes; sec Table 3 for components. 



at [Nem]/[Neii]~ 1, fall away rapidly. Among the 10 
galaxies in the sample with the lowest limiting band 
ratios L(7.7//m)/L(11.3/xm), all show evidence of low- 
luminosity LINER or Seyfert activity. In these systems, 
the relative strength of the short wavelength PAH emis- 
sion bands is suppressed by up to a factor of 10, com- 
pared to star formation dominated systems. An example 
of such a spectrum is shown in Fig. 15, for the LINER 
galaxy NGC 1316, with the predicted spectrum for the 
sample's median band ratios overlayed. 

One effect which could lead to the observed behavior is 
ionization state. Since the 11.3/im feature is thought to 
be produced by neutral PAHs, whereas the 7.7 /zm fea- 



ture arises primarily from PAH cations (e.g. Allaman- 
dola et al. 1999; Li & Draine 2001), a change in the 
neutral fraction could impact their relative strengths. 
Harder radiation fields are also more efficient at destroy- 
ing PAH grains, an effect which may dominate the ob- 
served trend, if the grain or molecular carriers of the 
7.7 fim and 11.3 /zm bands are different, and undergo 
photo-destruction via ionization or dissociation at sub- 
stantially differing rates. The 11.3 /im feature can be pro- 
duced by single-photon heating of PAHs that are larger 
than those that are effective at producing 7.7 /im emis- 
sion (e.g. Schutte et al. 1993; Draine & Li 2006). The 
observed decrease in L(7.7 /im) / L ( 1 1 . 3 /xm) may then be 
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Fig. 14. — The variation of inter-band strength ratios of the PAH features at 7.7/mi relative to the 11.3 /an band with the line ratio 
[Nciii]/[Nc n], an indicator of the hardness of the radiation field. Triangles (red) indicate galaxies with AGN (Seyfert or LINER) types, 
whereas squares (blue) are galaxies with H n-dominated nuclei. Upper limits are 2a. The dotted line shows the median of the H Il-dominated 
sources. Four sources (two AGN — magenta triangles — and two H II nuclei — cyan squares, labeled) are plotted with additional values 



obtained from spectra extracted over small 9 
values. 
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Fig. 15. — An example of the peculiar low-luminosity AGN 
PAH emission spectrum in barred LINER, galaxy NGC 1316, with 
the full decomposition shown as the solid green line. The PAH 
spectrum predicted from the strength of the 11.3 /imfeatures, and 
adapted to the SINGS median band ratios, is shown in blue, with 
a continuum-subtracted version below. 



the result of selective destruction of those PAHs small 
enough to emit at 7.7 /im. Low-metallicity dwarfs with 
the highest [Ne in] /[Neil] ratios appear to have roughly 
normal band ratios (when PAHs are detected), suggest- 
ing that AGN must play a role. 

Further evidence comes from a related trend in spa- 
tially resolved spectral maps of large SINGS spiral galax- 



ies with AGN nuclei. The 7.7/im/11.3/im band ratio 
in many cases decreases noticeably in the nucleus, com- 
pared to larger apertures which sample star-forming re- 
gions in the inner disk. An example of this effect for 
four galaxies chosen at d < 8 Mpc with both bright in- 
ner disk and nuclear emission is illustrated in Fig. 14. 
Ratios found in the default apertures are connected 
by solid lines to those recovered in small 9.25" x 9.25" 
(roughly 300 pc) apertures centered on the nucleus. For 
the two systems with AGN (NGC 5194 and NGC 7331), 
the short wavelength PAH features show further suppres- 
sion in their nuclear spectra, whereas the Hil galaxies 
(NGC 4631 and NGC 5713) move to slightly higher values 
of X(7.7/itm)/L(11.3/im). Such spatially-resolved varia- 
tions will be considered in more detail in a subsequent 
paper, but are consistent with the interpretation that the 
hard radiation fields near AGN have a strong impact on 
the strength of the short wavelength PAH bands. 

Kaneda et al. (2005) note an unusually weak 7.7/zm 
feature in IRS spectra of three X-ray-bright elliptical 
galaxies, as well as SINGS elliptical NGC 4125 (which 
has the fourth lowest PAH band ratio in Fig. 14). They 
suggest that unusual ISM properties in evolved ellip- 
tical galaxies could be responsible. However, two of 
these three ellipticals have established optical AGN ac- 
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tivity (NGC 2974, a Seyfert 2, and NGC 4589, a LINER), 
and our own recently obtained nuclear optical spectra 
of the third, NGC 3962, indicates that it also hosts an 
unambiguous AGN. This suggests that the AGN en- 
vironment, rather than the global morphological type, 
leads to the unusual emission characteristics. Elliptical 
galaxies, with little contamination from the PAH emis- 
sion associated with ongoing star formation, may sim- 
ply offer the most favorable systems for detecting this 
unusual spectrum. Consistent with this, Brandl et al. 
(2006) found a large scatter, but no significant variation 
of L(7.7 /xm)/L(11.3 //m) over a wide range of excitation 
in starburst nuclei. 

An alternative explanation is that the AGN is unre- 
lated to or only partially responsible for the change in 
the PAH spectrum. This would be the case if AGN clas- 
sification acts merely as a surrogate of low star forma- 
tion intensity, since strong nuclear star formation can 
overwhelm the AGN diagnostic lines. Cesarsky et al. 
(1998) found a depressed 7/xm band in differential ISO- 
CAM spectra within M31, suggesting this could be due 
to the weak UV field. In Spitzer samples of ^50 passively 
evolving early- type galaxies (Bregman et al. 2006; Bres- 
san et al. 2006), only two quiescent sources show PAH 
emission, one of which, NGC 4550, hosts a bright LINER 
nucleus and exhibits the same strong 11.3 /xm feature we 
report. The other, NGC 4697, is without an apparent 
AGN, yet also shows a moderately depressed 7 /xm PAH, 
though Bregman et al. conclude that after removing sil- 
icate emission, a spectrum similar to the diffuse ISM of 
the Galaxy is recovered. Future studies will better con- 
strain the role of AGN and weak star formation intensity 
in the modification of the PAH emission spectrum. 

If AGN dominate the effect, the absence of PAH emis- 
sion in nearly all quiescent galaxies suggests that, es- 
pecially in systems with little ongoing star formation, 
the AGN itself could excite PAH emission, supplying 
UV radiation dominated by Lyman-a and H2 Lyman- 
Werner bands in the extended X-Ray dominated regions 
(Meijerink & Spaans 2005), with a hardness spectrum 
quite unlike that typically found in environments har- 
boring young stellar populations. In addition, an AGN 
may emit X-rays sufficient to destroy small PAHs di- 
rectly even at distances of kpc (Voit 1992). Because 
larger PAHs are expected to be less susceptible to de- 
struction by X-rays, the presence of an AGN can alter 
the PAH emission spectrum, even if the PAH excita- 
tion itself is dominated by starlight from an older stellar 
population. Clearly, the absolute luminosity density of 
the AGN source is also an important factor regulating 
the PAH spectrum, since more luminous AGN either de- 
stroy all PAH grains (e.g. Voit 1992; Lutz et al. 1998), 
or produce mid-infrared continuum photons in such co- 
pious quantities that the PAH bands are diluted away 
(Siebenmorgen et al. 2005; Hao et al. 2005; Sturm et al. 
2005). This effect was exploited by Genzel et al. (1998) 
to differentiate between AGN and starbursts at high lu- 
minosity. Although grain size modification may operate 
in many contexts, the unusual band ratios observed in 
low-luminosity AGN nuclei (Z(7.7/xm)/Z(11.3/xm) < 2) 
could offer an efficient way to detect such weak active 
nuclei, when the UV field is not dominated by star for- 
mation, and dust absorption hinders optical classifica- 
tion. Though it is not yet known over what regimes this 
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Fig. 16. — The variation of the 17 /tm/11.3 (im PAH lumi- 
nosity ratio with oxygen abundance. Red triangles are sources 
with LINER/Seyfert classifications, blue squares have Hn classifi- 
cations. 



effect operates, the possibility that AGN sources could 
directly excite PAH emission suggests considerable cau- 
tion when using the PAH bands as direct indicators of 
star formation in systems known to harbor active nuclei. 

Care must be taken to quantify the band ratios which 
constitute an unusual PAH spectrum of this type, since 
the 11.3/xm feature, while not the most luminous, has 
the highest specific luminosity per wavelength interval, 
and does not typically suffer continuum dilution from 
starlight as the shorter wavelength PAH features can. 
At low signal and/or with strong stellar continuum, the 
11.3 /mi feature can remain the most easily detectable 
PAH band, and appear to dominate despite normal band 
ratios. An example of this phenomenon is found in the 
SINGS dwarf galaxy NGC 1705 (Cannon et al. 2006). 

5.6.2. PAHs and Metallicity 

That the absolute strength of PAH emission can de- 
pend on metallicity has recently been well established 
(Engelbracht et al. 2005; Madden et al. 2006; Wu et al. 
2006; O'Halloran et al. 2006; Rosenberg et al. 2006), with 
indications of a threshold 12+log(0/H) ~ 8.1, or approx- 
imately Zq/A, below which PAH emission is suppressed. 
From detailed SED fitting of the SINGS sample, Draine 
et al. (2006, in prep) find that the fraction of the dust 
mass contributed by PAHs is suppressed below a simi- 
lar threshold. Most studies rely on Spitzer broadband 
photometry (e.g. 8/xm/24/xm) to assess absolute PAH 
strengths. Although broadband 8 /im flux does correlate 
well with PAH strength (see §5.5), the range of mid- 
infrared continuum properties can overwhelm even sig- 
nificant changes in the relative PAH luminosity — see, 
e.g., the large variation in 7.7 /xm PAH to 24 /xm luminos- 
ity in Fig. 11. Broadband methods are also not sensitive 
to changes in the form of the PAH emission spectrum. 

Figs. 16 and 17 illustrate two examples of the behavior 
of PAH emission with changes in nuclear oxygen abun- 
dance: the relative strengths of the 17/xm and 11.3 /xm 
PAH bands, and the total integrated PAH luminosity rel- 
ative to the total infrared. Oxygen abundances for the 
SINGS galaxies have been determined by Moustakas et 
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Fig. 17. — The variation of the total PAH luminosity relative 
to the total infrared (3-1100 /im) luminosity. The threshold abun- 
dance 12+log(0/H) = 8.1 is indicated. Symbols are as in Fig. 16. 



al. (2006, in prep.), based on a combination of nuclear 
and circumnuclear optical spectroscopy and observations 
of individual Hn regions culled from the literature, and 
placed on the common 12+log(0/H) strong- line calibra- 
tion scale of Pilyugin & Thuan (2005). 

The central spectra of SINGS galaxies probe much 
higher oxygen abundance than the \ow-Z sample of En- 
gelbracht et al. (2005), yet, even at these relatively high 
metallicities, up to ~ 1.5Zq, there is a definite trend 
in the inter-band strength ratios, with higher metallicity 
leading to increased L(17/zm)/L(ll/im). Though low 
metallicity systems do not tend to harbor AGN (Kauff- 
mann et al. 2003), the presence of a weak AGN in the 
nucleus clearly has a strong impact on the behavior of 
L(17/J,m)/L(ll /jm), with the highest ratios observed 
only in systems with active nuclei. 

Since the 15-20 /jm PAH bands are hosted by larger 
grains than the shorter wavelength features, the trend 
with metallicity is suggestive of a scenario in which for- 
mation of the large grains that contribute the 17 ^m band 
is enhanced in the presence of higher metal abundance. 
This effect could act in combination with the destruc- 
tion of smaller PAH grains in AGN environments, since 
it is likely that the mechanism that suppresses the short- 
est wavelength PAH features in AGN has a similar, if 
smaller effect on the 11.3 /im feature. 

Low metallicity systems also have harder radiation 
fields due to a combination of two effects: 1) hotter 
massive stars form in low metallicity environments (an 
effect dominated by the metallicity sensitivity of line- 
driven stellar winds, Martin-Hernandez et al. 2002), and 
2) dust production is inhibited, leading to reduced shield- 
ing by dust grains. This could also impact the PAH 
band ratios of Fig. 16. However, since both the 11.3 /tin 
and 17 /mi band are contributed by neutral PAHs, the 
ionization state cannot drive the trend in their ratios. 
Likewise, selective destruction of the smaller grains in 
the harder fields would result in the opposite trend - 
higher L(17/mi)/Z/(11.3/im) at lower metallicity. This 
leaves the inhibited formation of the largest grains at low 
metallicity as the most likely dominant effect. Clearly, 
however, metallicity variation alone cannot explain the 



full range of L(17/zm)/L(11.3/i,m) in Fig. 16, since at 
Z>Zq/2, AGN consistently have the largest ratio. 

The behavior of the total PAH luminosity, integrated 
over all bands from 5-19 /im, relative to the total infrared 
luminosity in Fig. 17 shows a similar trend, with high 
fractional PAH contribution to the total infrared (> 10%) 
occuring only at metallicities above 12+log(0/H) ~ 8.1. 
There is also a marked separation between galaxies with 
weak AGN and H n classifications. At a given metallic- 
ity (even super-solar), AGN sources emit a much smaller 
fraction of their central infrared luminosity as PAH emis- 
sion. Coupled with the anomalous PAH band ratios seen 
in the most extreme such sources (Fig. 14), a picture of 
weak AGN profoundly influencing the absolute strength 
and shape of the PAH emission spectrum emerges. 

6. SED MODELS IN THE HIGH-REDSHIFT CONTEXT 

Rather than directly modeling the physical processes 
and resulting mid-infrared emission properties of star- 
forming galaxies, most broadband infrared studies of 
moderate to high redshift galaxies have utilized infrared 
templates, constructed from small samples of local galaxy 
spectra. Since the PAH bands are so important energet- 
ically, and impart the 3-19 /jm spectrum of star-forming 
galaxies with features that contrast strongly with each 
other and with the underlying dust and stellar contin- 
uum, they are a vital component of any infrared SED 
template. A small number of bright, very well stud- 
ied systems with high-quality ISO spectra (e.g. Arp 220, 
M82, NGC253, Mrk231; Rigopoulou et al. 1999; Sturm 
et al. 1996, 2000) have been used extensively (e.g. Elbaz 
et al. 2002; Webb et al. 2006; Reddy et al. 2006). The av- 
erage ISOPhot spectrum of Lu et al. (2003) has also seen 
wide used, having been integrated into the SED library 
of Dale et al. (2001) and others. Even physically realistic 
SED models (e.g. Silva et al. 1998) are often used with 
the mid-infrared substituted by scaled template spectra. 

The importance of quantifying not just the average 
shape of the mid-infrared spectrum of star-forming galax- 
ies, but the range of variation in that shape, is well il- 
lustrated by the source count models of Lagache et al. 
(2004). To model deep MIPS 24 /mi source counts (Pa- 
povich et al. 2004), Lagache et al. applying a modest 
tilt to the ISO PAH template spectrum to enhance the 
power at 15 /im. With this modified SED, the source 
counts, which could not be reproduced by pre-existing 
models, came into very good agreement. This pertur- 
bation played a large role in the interpretation of the 
redshift distribution and star formation history implied 
by the broadband survey sample. Such sensitivity to 
the detailed form and shape of the mid-infrared emission 
spectrum underlines the urgent need to quantify vari- 
ations in this spectrum, and establish the role of such 
natural variations in the interpretation of resolved and 
unresolved survey results. 

In Fig. 18, five noiseless model templates from the spec- 
tral decomposition technique of § 4 are shown — an av- 
erage spectrum weighted by the global total infrared lu- 
minosity, and four sub-sample averages drawn from the 
empirically-defined regions in Fig. 13, formed by lines 
in the L(6.2/tm)/L(7.7/tm) vs. L(11.3/im)/L(17/im) 
space. These boundaries were chosen to ensure ade- 
quate sampling in each region, and to bracket a modest 
range of variation within the sample. Only galaxies with 
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Fig. 18. — Noise- free MIR, templates, produced by a TIR 
luminosity- weighted average of all spectra (black), and from sub- 
sample means for each region in Fig. 13 (color). The units are arbi- 
trary vl v , with a constant offset of 2.5 units between spectra. Mean 
matched 70 fim and 160 /an intensities are also shown. The lower 
panel contains the normalized residuals (model-spectrum) /model, 
averaged over the sample, with 20% trimmed from the extremes of 
the distribution. 



reliable detection in all four main PAH bands are in- 
cluded (eliminating all sources dominated by stellar con- 
tinuum). Emission lines are not removed. The region- 
averaged spectra corresponding to the lower two quad- 
rants in Fig. 13 have the largest contribution from weak 
AGN sources. The normalized residuals are also shown, 
and are typically only a few percent, with some struc- 
ture near known artifacts (e.g. fringing at 20-23 /im). 
The luminosity- weighted average has the steepest 30 /im 
continuum, an indication that the most luminous sources 
have the warmest thermal SEDs. The 70 /im and 160 /im 
matched intensities (as described in § 3.2) are shown for 
each template spectrum. 

The full sample was used to investigate the redshift 
response of the Spitzer bands, defined as 



R(z) 



J 9xh{z) dX 

r>20f-L7n 



SsTh(z)dX 



(5) 



where g\ is the effective filter transmission function for 
the three broadband Spitzer imaging filters of interest: 
IRS Peak-Up Blue (16 fun) and Red (22 /mi), and MIPS 
24 /im. R{z) quantifies the expected broadband response 
to a galaxy of a given spectral appearance. The nor- 
malization limits were chosen to ensure that the PAH 
spectrum does not suffer artificial variations due to the 
large range of slopes in the 20-30 /im dust continuum. 

In Fig. 19, R(z) is plotted relative to its sample mean 
at each redshift, indicating a wide range in broad pho- 
tometric response arising from normal variations in the 
PAH spectrum, and, to a lesser degree, the underlying 
continuum. In the highest luminosity quartile, the range 
of variation is reduced by roughly 30%, and there is a 
waist of diminished variation (a factor of 1.4) at z ~ 0.6, 
in the relatively stable rest-frame continuum between the 
11 fim and 17 /im bands. At the low and high redshifts 
ends, variations in the continuum shape dominate the 
response. The median range in the 10%— 90% trimmed 
sample is roughly a factor of two (1.6 in the most lumi- 
nous quartile). At redshifts where the region bracketing 
the 7.7 [im and 11.3 /im features shifts through the bands, 
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Fig. 19. — The redshift response of three broadband Spitzer 
filters (IRS PU Blue at 16 /im, IRS PU Red at 22 (im, and MIPS 
24 fim) to the range of sample spectra is shown, relative to the 
mean response at each redshift. The global luminosity-weighted 
mean relative response is shown as the dashed line (the mean PAH 
template spectrum would produce the flat line R(z) = 1.0). The 
light gray region encloses the 10%— 90% range of variation within 
the full sample, and the dark gray region encloses a similar 10%- 
90% range for the most luminous quartile of the sample (£tIR > 
2.6 X 10 10 Lq). The redshifts where the peak wavelengths of the 
main PAH features are shifted to the 50% power and transmission- 
weighted central wavelengths of each broadband filter are indicated 
by shaded regions for the main PAH bands at (from left to right, 
red to blue): 17 fim (dotted, red), 11.3 fan (long dash, red-blue), 
7.7 /im (dash-dot, blue-red), and 6.2 /tm (dash-dot-dot-dot, blue). 



the range of variation increases significantly. 

At luminosities above Ltir ~ 10 n L©, other MIR emis- 
sion components can begin to dominate over PAHs. In 
particular, the hot dust emission of embedded AGN can 
contribute the bulk of the mid-infrared power and dilute 
away any PAH emission present, and very strong sili- 
cate absorption can have a profound effect on the shape 
of the MIR continuum. A fundamental question in the 
modeling of deep infrared survey results is what com- 
bination of features constitute a useful suite of galaxy 
template spectra. There is no guarantee that distant 
galaxies forming in the early universe will exhibit the 
same infrared properties as local galaxies, even if care- 
fully sub-divided by total infrared luminosity, metallic- 
ity, or other systemic physical properties. Because PAH 
emission is detected or inferred at ever higher redshifts in 
systems with a wide variety of infrared properties, apply- 
ing the complete range of locally observed PAH behavior 
is the only option which fully brackets the known uncer- 
tainties which result from natural variations in the PAH 
emission spectrum. 

7. SUMMARY AND CONCLUSIONS 

We have presented a new suite of 5-38 /im spectra from 
the inner ~l/2 arcmin 2 (roughly 1-10 kpc 2 ) of a broadly 
selected sample of nearby, predominantly star-forming 
galaxies, many of which host weak AGN. The spectra are 
dominated in most cases by PAH emission bands limited 
to the range 5-19 /im, stellar photospheric emission of 
varying strength, smooth thermal dust continuum, and 
occasional absorption by silicate grains. 

The low-resolution spectra of star-forming galaxies can 
be successfully reproduced using a simple physically- 
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motivated decomposition model with a small number 
of components, including starlight, thermal dust contin- 
uum, and line and dust emission features, all absorbed 
by single parameter dust extinction. The decomposition 
method is also successful at recovering unresolved emis- 
sion lines blended with dust features. 

Weak dust emission features at 6.69 /zm at 18.92 /im 
are reported for the first time in cxtragalactic sources, 
and many of the weak features formerly seen only in ISO 
spectra of one or two of the brightest starburst galax- 
ies are found commonly among our sample. The 17/im 
complex correlates strongly with other PAH bands, sup- 
porting the identification of this new feature with PAH. 

Compared to spline-based or other simple techniques 
for estimating the strengths and continuum underlying 
PAH features, the decomposition method recovers signif- 
icantly more of the extended feature power. Spline meth- 
ods agree reasonably well within a simple scalar factor 
for the two main bands in clean continuum regions at 
6.2 /im and 11.3 /im, but recover much less of the total 
feature power, and with significantly more scatter, for 
the blended complexes at 7.7 /zm and 8.6 /zm. 

We have formulated a new extinction profile for use in 
the 5-38 /im region. Measurable dust absorption at mid- 
infrared wavelengths by silicate grains is found to occur 
in ~l/8 of the sample targets. Most normal star-forming 
galaxies appear to suffer less than Ay ~ 3 magnitudes of 
foreground screen silicate extinction, averaged over their 
central few square kiloparsecs. 

Matched Spitzer photometry at 24, 70, and 160 /im, 
tracing the bolometric dust emission, highlights the re- 
markable power emitted solely in the PAH bands, which 
in aggregate can contribute up to 20% of the total 3- 
1100 /im infrared luminosity, with a typical value of 10%. 
The PAH features vary by a factor of 10 relative to 24 /tm, 
with indications that both stochastically and thermally 
heated grains contribute to the continuum there. 

The individual PAH bands show considerable variabil- 
ity — factors of 2-5 in internal strength ratios in the 
sample — with some features varying more than others. 
The 17/im band in particular exhibits variation which 
tracks with metallicity, such that systems with higher 
oxygen abundance are more likely to have stronger 17 /im 
emission relative to the shorter features. Since the longer 
wavelength PAH emission is expected to arise from larger 
PAH grains, this could indicate that the size distribution 
of the PAHs depends on the metal content of the ISM, 
even at relatively high metallicities reaching 1.5Z Q . 

Contrary to the effects of luminous AGN, which com- 
pletely destroy PAH grains or mask their emission bands 
with hot dust, the presence of low-luminosity AGN in 
the nucleus of a galaxy has a profound effect on the form 
of its PAH spectrum, leading in some cases to peculiar 
emission with very weak or absent 6.2 /im, 7.7/zm and 
8.6 /zm bands, suppressed in strength by up to a factor 
of 10. Only sample galaxies with low-luminosity LINER 
or Seyfert classifications from their optical nuclear spec- 



Allamandola, L. J., Hudgins, D. M., & Sandford, S. A. 1999, ApJ, 
511, L115 

Allamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. 1985, 

ApJ, 290, L25 
Armus, L., et al. 2004, ApJS, 154, 178 

Barvainis, R., Antonucci, R., & Helou, G. 1999, AJ, 118, 645 



tra produce these exceedingly unusual band ratios, rais- 
ing the possibility that the weak AGN itself can directly 
modify the PAH grain size distribution, and even serve 
as the excitation source for PAH emission. If so, the ab- 
normal PAH strengths could be used to detect such weak 
AGN systems in dusty galaxies where optical diagnostics 
are unavailable or unreliable. Alternatively, low star for- 
mation intensity in the centers of many AGN may play 
a role in modifying the form of the PAH emission. 

Both higher metallicity and the presence of an AGN 
shift power from the short to long wavelength PAH 
bands, most likely by rebalancing the distribution of 
grain sizes by favoring formation of larger grains, and 
selectively destroying the smaller grains which emit the 
bulk of the shorter wavelength PAH features. The in- 
tegrated luminosity of the PAH bands, relative to the 
total infrared, also shows an overall increase with increas- 
ing oxygen abundance, but this trend does not hold for 
galaxies hosting AGN, which emit a much smaller frac- 
tion of their total infrared luminosity in the PAH bands 
(a few percent), even at relatively high metallicity. Both 
the unusual PAH spectrum, and the diminished PAH lu- 
minosity in these sources caution against the direct use 
of the absolute PAH band strength as an indicator of star 
formation rate in systems known to harbor AGN. 

Taking into account a trimmed range of variations in 
the PAH band ratios and the shape of the mid-infrared 
continuum seen in the sample, we find that, compared to 
the common assumption of an invariant PAH emission 
template used to model high-redshift sources, at least 
a factor of two variation in broadband photometric re- 
sponse is imparted at 16 /im, 22 /tm, and 24 /im, at all 
rcdshifts z <2.5. At higher total infrared luminosities, 
above 2.6 x 10 10 L Q , the effective range is smaller by 
roughly 30%, but still substantial. While it is likely that 
as galaxy luminosities approach 10 12 io, other effects, 
such as the growing contribution of obscured AGN to the 
typical infrared SED, will dominate the uncertainties in 
modeling deep mid- infrared source populations, utilizing 
a single, fixed PAH template will impart factor of two 
or larger uncertainties due solely to natural variations in 
the form of the mid-infrared PAH emission spectrum. 
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